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Abstract 

In this letter, a theoretical method for the analysis of diffusive flux/current to limited scale self- 
affine random fractals is presented and compared with experimentally measured electrochemical 
current for such roughness. The theory explains the several experimental findings of the temporal 
scale invariance as well as deviation from this of current transients in terms of three dominant 
fractal parameters for the limited- length scales of roughness. This theoretical method is based 
on limited scale power-law characterization of the interfacial roughness spectrum and the solution 
of diffusion equation under the diffusion-limited boundary conditions on rough interfaces. More 
broadly, these calculations challenges the earlier belief that the anomalous behavior is solely de- 
pendent on fractal dimension of roughness and highlight the potential to apply this equation for 
the scale invariant roughness determination. Finally, the validity of theoretical result is tested with 
extensive experimental data. 

PACS numbers: 66.10.Cb, 82.45.Yz, 47.27eb 



Realistic surface roughness has hmited-lenerth scales of irregularities and are frequently 

n p nti 

characterized as self-afiine fractals pj, |2, 13|, 1^- Diffusion-limited processes on such in- 
terfaces show anomalous behavior of the reaction flux. Some of the diverse realizations 
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131 and electrochemistry 




of diffusion-limited processes in physicalphenomena are: Spin relaxation [sl, fluorescence 
quenching , heterogeneous catalysis[7|, |8j, e nzy me 
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The diffusion of the reactant from a bulk medium towards an interface where the reactants 



either loose their activities or are trans: 
areas of science 
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product is a common problem in diverse 
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2J]. These 



processes are experimentally approximated as power-law relation of reaction rates/flux, J{t), 
in time {t) and is represented by the following relation: 



J(t) ~ t 
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where the exponent, /? depends on interfacial roughness. Theoretical justification for Ec 
was provided by De Gennes scaling result with /3 = {Dh — l)/2j^ and generalized form 
DeGennes analysed it for the problem of diffusion-controlled nuclear magnetic relaxation in 
porous media with fractal interfacial dimension Dh- Later, similar results were discussed 
for other diffusion-controlled situations, such as adsorption on a porous fractal catalyst 



in the context of flow of energy and mass through a fractal interface 
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for rough fractal 



15| and similar result 



electrode/electrolyte interfacial current under potentiostatic conditions 
for heat diffusion from a self-affine fractal boundary[ll|. Due to its simplicity and the lack 
of better alternative of Eq. ([1]) captured lots of attention and is used indiscriminately for 
interpreting large quantity of data 20|, l2l|]. These include cases where the range of roughness 



is too small to be taken as idealized fractals. It is commonly known that Eq. ([T]) is unable to 
include complete characterization of realistic fractal roughness and also does not reproduce 
subtle aspect of experimental data. 

I, present an analytical model that allows one to probe the anomalous time behavior 
of diffusive flux of reacting interface with roughness and to evaluate the consequences of 
limited length scales of roughness on the anomalous behavior. This theory explains the 
several experimental findings 20|, l2l|] of the temporal scale invariance of flux/current in 
terms of three dominant fractal parameters for the limited-length scales of roughness. 

The information about surface roughness enters in my theory through power spectrum 



2 



of roughness. The power spectrum of a reahstic surface (also called " approximate self-afEne 
fractal" ) is described in term of limited scales of wave- numbers {K) power- law function 
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I.e., 



C{K) ) = T\K\'^ " , for 1/L < \K\ < l/£. There are four physical parameters 



ractalitYj respectively; and T is the 



3|, [261,1271], its units are cm^^^-^ and 



of roughness in this framework, namely D^, i , L and T. is the fractal dimension, a 
global property which describe scale invariance property of the roughness- an anomalous 
behavior in flux and its time exponent is usually assumed to be function of this parameter; 
£ and L are lower and upper cutoff length scales of 
strength of fractal and related to topothesy of fractals 
T — > implies no roughness. The lower roughness scale is the length above which surface 
show fractal behavior. Also, such power-law spectrum may represents band-limited form 
of Ausloos-Berman's generalization [28] for multivariate Weierstrass-Mendelbrot function for 
an isotropic (statistically) rough surface. The moments of power-spectrum are related to 
various morphological features of rough surface viz. rms width(y'mo), rms gradient (-Y/m^), 
rms curvature (ym4) etc. The general moments of power spectrum (i.e., 2k-th moments, m2fc) 
are easily obtained for above mentioned power-spectrum and are important morphological 
characteristic of surface roughness. The general formula is: m2k = 
where 6k = S + k and 6 = — 5/2. 

Present theory of anomalous diffusive reaction rates on realistic self-afiine fractals is 
obtained using our general formalism jirl-fisl which show the diffusive fiux/current at random 
rough surfaces can be described in term of its power spectral density. Main approximation 
involved in derivation of general formalism is the truncation of solution at second order 
in ..face .oughnes. p.ofi.H Q. The total (averaged) flux/ou„.eM at the stationary. 
Gaussian random surface is given byfl?!. [l8|: 
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where D is the bulk molecular diffusion constant, Cg is the difference between the surface 



and bulk concentration and Aq is the area of surface around which rough surface fluctuates. 
The diffusion controlled reaction rates are related to potentiostatic current transients of an 
electrode undergoing fast charge transfer. The reaction flux, J(t), is related to electrode 
current (/(t)) as: I{t) = ~nFJ{t), where n is the number of electron transfer in redox 
reaction {OsoiuUon + ne~ ^ RsoiuUon) and F is Faraday constant. 

The exact solution for the dynamic diffusive flux on an approximate self-affine surface 



(substituting above mentioned band- limited power law spectrum in Eq. ([2])) under diffusion- 
limited condition is obtained as 

_ DA,Cs f^^T f i-^' - L-^' r(^, Dt/e, Dt/L') \\ 
^^"v^V SDt + {Dty^s )) 

where 6 = — 5/2, V (a, xq, Xi) = T (a, xq) — T(a, xi) = 7 (a, Xi) — 7 (a, Xq), T{a, x,) 



and 7(a,Xj) are the incomplete Gamma functions [29|]. Equation ([3]) graphically analysed 
in Fig.l which show dependence of scaling region on i, Dh and T. The scaling region has 
very weak dependence on L so it is not shown in Fig. 1. 

Most of experimental data recorded are for the intermediate time regime i.e., L'^/D }^ t > 
£'^/D. The expression for intermediate time is obtained using time constraints in expansion 
of two incomplete gamma functions and retaining only leading orders. Final equation has a 
simple and elegant form as follows: 

^' V^t\ S7r\SDt (Dty*')) 

which capture the anomalous behavior of reaction flux. r(x) in Eq.Q is the gamma 
function [29!]. The small length scale of fractal roughness command the flux transient for 
roughness of limited scales, i.e., L/£ ~ 10. There is no dependence on the upper length 
scale(L) for fractality in Eq. (jlj). 

Equations ([3j) and (jlj) extends the conventional representation of the Cottrell current 
transient ( 1/Vi) on the planar electrode in electrochemistry {30I to the fractally rough 
electrode. This equation achieves a more realistic characterization of limited scale rough 
surface diffusive flux as it includes the fractal dimension dependent power-law as well as 
the contribution from the length scale and strength of fractality. The total mean flux is the 
summation of smooth surface flux and an anomalous excess flux due to fractal roughness 
or can also be looked upon as product of the l/y/i current and dynamic roughness factor 
(term inside parenthesis of Eqs. ([3]) and (jlj)). 

Nature of the plots also elucidates anomalous scaling behavior in the intermediate time 
regime. The current transients increased with decrease in small scale of roughness in the 
early time domain, which is similar but larger in the magnitude of slope as compared to the 
planar 1 / ^/t response. It follows a power law behavior in the intermediate time regime which 
merges with large time 1 / \/i behavior. As the range of roughness increases, the roughness 
factor and total current output also increases simultaneously on small time scales. So one 
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FIG. 1: Illustrates the effect of three dominant fractal parameters, i.e., Dh, i and T, on the 
anomalous scaling behavior of reaction rates. These theoretical results are plotted as graphs to 
show an approximate power-law dependence of the flux / current on time. These calculations break 
the earlier beliefs based on idealized fractal models that the exponent of anomalous diffusion region 
purely depends on fractal dimension of roughness, a, Anomalous region dependence on the fractal 
dimension of roughness. The solid line is generated using f = 0.6/xm, T = 1.3x10^^ {arb. units) and 
Dh = 2.45, 2.40, 2.35 and 2.30 from above. Ratio of L/i is kept fixed (=100) in all calculations, 
b, The effect of lower scale of roughness on the anomalous region is clearly depicted in the 
double logarithmic plot of current transients. The solid line is generated using Dh = 2.4, T = 
1.3x 10~^(ar6. units) and I = 0.2 fim, 0.4 fim, 0.6 iim and 0.8 nm from above, c. Shows the effect of 
strength of fractal on the anomalous region or slope of double logarithmic plot of current transients. 
The solid line is generated using Dh = 2.4, £ = 0.4 /xm and T = Tq = 1.3 x 10~^ {arb. units), 
1.5xro, 2xro and 2.5xro from above. Other fixed quantities in calculations are: macroscopic areas 
{Ao = 0.1 cm^), diffusion coefficient (D = 5 x 10~^ cm^/s) and concentration {Co = Cr = 5 mM). 
These figures clearly demonstrate that the exponent of the anomalous region is dependent on all 
three fractal roughness parameters. 



can say that the total current output at small times is dependent upon the lower cutoff 
length scale I. No such impact of £ has been observed at large time scale as this regime is 
controlled by the width of interface ttlq. The width of interface is a strong function of L and 
a weak function of £. Another important feature that can be verified from these graphs is 
that there is no such sharp outer cut-off time but the inner cut-off time decreases with the 
decrease in the lower length scale of roughness or with the increase in roughness factor. 

What is remarkable about equation (jl]) is that it governs the diffusive flux of the diverse 
set of roughness features. Its validity is not only tested, in this letter, for various range 
of roughness but also for the magnitude of roughness factor too. The conventional scaling 
equation (1) cannot explain deviation from linear behavior in log-log plot which is often seen 



in data 



2l| and is easily captured by this theory. Figure (2) has several current transient 



curves and experimental data 



21| on rough electrodes with limited length scales of fractality. 



Data are analyzed keeping in mind the experimental information about roughness factor that 
is available about their roughness. Two lengths I and L are approximately identified from 
STM image of surface 211]. Identification of Dh-, ^ and T is achieved by the minimization of 
deviation from experimental time-dependent data of diffusion-limited current. 

Surprisingly, limited order perturbation analysis is able to capture the features of large 
roughness too. Though one expects that the scaling results would be seen only in large 
roughness form of theory to match with assumptions of Eq. ([1]) [3] ■ Most important of all, 
this work shows an intermediate anomalous power-law form for time above inner transition 
time {ti). This suggests that this theory have extended range of validity much beyond 
expectations which is also seen in comparison with experimental results in Fig. ([2]). Another 
feature of this theory is that ti decrease with decrease in ^. Inner transition time is evaluated 
equating small time expansion 

and Eq. (jl]) at t = U. 

Results developed in Eqs.Q and (jl]) are based on statistical models for the current 
on random surface fractals. Advantage of such formulation is that it is based on four 
statistical parameters: Dh-, T, i and L. Most cases, we do not have independent information 
about fractal parameters as experimental studies rarely measure roughness power spectrum. 
In some cases, small and large length scales of roughness is characterised from SEM or 
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FIG. 2: Comparison of model predictions from Eq. ^ and current transients data of Ocon et al 21]. 
Electrodes in these diffusion-limited current experiments used were nano scale electrodispersed 
gold wire with surface layer of thin columnar gold 2l| . There roughness is changed by ageing of 



these columnar structure gold deposits by a slow reorganization processes. Ageing may lead to 
increase in £ and decrease in Dh of rough surface. Rough surface of a wire electrode is imagined 
to be randomly fluctuating surface around a macroscopic plane as the curvature contribution of 
macroscopic wire geometry is insignificant, a, The solid lines are generated using: (1) Dh = 2.42, 
£ = 11 nm, T = 2.0 x 10^'^{arb. units), Aq = 0.93 cm?, and experimental data points (A). (2) 
Dh = 2.19, £ = 400 nm, T = 4.0 x 10~^ {arb. units), Aq = 0.073 cni^ and experimental data 
points (V). (3) Dh = 2.02, £ = 750 nm, T = 5.0 x 10~^{arb. units), Aq = 0.068 cm^ and 
experimental data points ((>)• b. The solid line is generated using: (1) Dh = 2.52, ^ = 20 nm, 
T = 1.94 X 10^'^ (ar6. units), Aq = 0.25 cm^ and experimental data (a). (2) Dh = 2.33, ^ = 86 nm, 
T = 1.04 X lQ~^{arb. units), Aq = 0.085 cm^ and experimental data (t). (3) Dh = 2.32, 
£ = 87 nm, T = 1.32 x 10^^(ar6. units), Aq = 0.060 cm^ and experimental points (♦). 

STM image of surface. Similarly, the knowledge of surface roughness factor or width of 
roughness or both can help to fix these statistical parameters. The remaining unknown 
parameters can be obtained by minimizing variance of experimental data for current from 
theoretical values. Knowledge of these four parameters one can predict various roughness 
feature of roughness profile like: mean roughness factor (-R*), root mean square (rms) width 
of roughness {h or ^rrio)), rms gradient {^Jm^) and inverse rms curvature {r^ = l/y/m^). 



The roughness factor is a function of mean square gradient m2, i.e., R* ~ A/7rm2/2 for 
large roughness, which in turn is a strong function of i and a weak function of L. Figure 
2 compares our theory for current transient with rough gold deposit on wires. Equation 
(jl]) is also obeyed by other experiments on surfaces like replica or gold masking of surfaces 
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like fractured steel, dental surface and liquid-liquid interface. This theory predicts various 
roughness and marpho logical features and tabulated in Table 1 along with tj. 



Table 1 Predicted morphological parameters & ti 



Data 






i?* h (/im) 




Tc (nm) 


U (ms) 


Fig 2 a (1) 


2.4 


2.42 


50 


1.79 


39.90 


0.40 


0.023 


Fig 2 a (2) 


2.2 


2.19 


3 


2.19 


2.25 


269.3 


0.350 


Fig 2 a (3) 


2.0 


2.02 


1.1 


0.50 


0.34 


3344.8 


0.743 


Fig 2 b (1) 


2.52 


2.52 


100 


4.87 


79.97 


0.35 


0.143 


Fig 2 b (2) 


2.24 


2.33 


20 


4.65 


15.96 


7.99 


0.522 


Fig 2 b (3) 


2.32 


2.32 


20 


4.76 


15.95 


8.12 


0.546 



Dh^ and Dh^ are calculated using Eci^(l) and ([4]), respectively. Roughness factor i?* experimentally mea- 
sured from voltammetric experiment [21] and is used to constraint value of T parameter in our calculations, 
Predicted width of interface(/i), root mean square gradient (^7712) and average radius of curvature(rc) and 
inner crossover time(ii). Upper cutoff length (L) is kept constant in all calculations i.e., L = 3 fim. 

The central result, equation (jlj), constitutes an elegant and simple test that any data must 
pass to be called limited scale self-affine fractals. This model gives a very good description 
of large quantity of data which have the scaling region as well as region which deviate from 
it. It is important to note that the slope of scaling region does not purely depend on fractal 
dimension alone but also on lower cutoff length scale of roughness and strength of fractality 
too. Another important observation is that the intermediate time expression is excellent for 
the current transient in intermediate and long time regimes (i.e., t > ti). 

This formulation opens new avenues towards an understanding of diffusion -limited re- 
action rates at realistic rough surfaces with limited scales of fractality and inverse problem 
of obtaining fractal parameters. We demonstrate that the fractal dimension is not in it- 
self sufficient to understand this problem. This work unravels the connection between the 
anomalous intermediate power-law regime exponent and the geometric parameters of lim- 
ited scales of fractality i.e. Dh, T. I believe that this will provide insights into many 
other diffusion-limited processes on random interfaces and the generality of this approach 
will provide logical extension to several complex situations, like various transient techniques 
in electrochemistry, finite reaction rates at interface, coupled homogeneous reactions with 
heterogeneous interfacial reaction etc. As more data becomes available with clear power 
spectral characterization of roughness, quantitative analysis should provide an improved 
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way to understand roughness of several surfaces. 
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